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Abstract. We consider a generalization of the spin-boson model in which two 
different two-level systems are coupled to an oscillator, under conditions where the 
oscillator energy is much less than the two-level system energies, and where the 
oscillator is highly excited. We find that the two-level system transition energy is 
shifted, producing a Bloch-Siegert shift in each two-level system similar to what would 
be obtained if the other were absent. At resonances associated with energy exchange 
between a two-level system and the oscillator, the level splitting is about the same as 
would be obtained in the spin-boson model at a Bloch-Siegert resonance. However, 
there occur resonances associated with the transfer of excitation between one two-level 
system and the other, an effect not present in the spin-boson model. We use a unitary 
transformation leading to a rotated system in which terms responsible for the shift and 
splittings can be identified. The level splittings at the anticrossings associated with 
both energy exchange and excitation transfer resonances are accounted for with simple 
two-state models and degenerate perturbation theory using operators that appear in 
the rotated Hamiltonian. 



PACS numbers: 32.80.Bx,32.60.+i,32.80.Rm,32.80.Wr 
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1. Introduction 

The coupled quantum system consisting of a two-level system interacting with a 
harmonic oscillator provides a model that has been focus of a large number of studies 
over recent decades. A reduced (Rabi Hamiltonian) version of the problem was published 
by Bloch and Siegert in studies of the interaction of a dynamical magnetic field 
perpendicular to a static magnetic field with a spin system pQ, and was later applied to 
the problem of atoms interacting with an electromagnetic mode [2] . The more complete 
model (spin-boson Hamiltonian) was introduced much later by Cohen- Tannoudji and 
collaborators [3]. 

At the most basic level, the interaction of the two-level system with the oscillator 
produces a shift in the two-level energy (the Bloch-Siegert shift). Energy exchange 
between the two-level and oscillator is allowed when the shifted two-level system energy 
becomes resonant with an odd number of oscillator quanta, which results in level 
anticrossings at these resonances (Bloch-Siegert resonances). If the characteristic energy 
of the oscillator hu is much less than the transition energy AE, then many oscillator 
quanta are required to match the shifted (dressed) transition energy. In such a limit, the 
model is considered to be in the multiphoton regime; which is a topic of current interest 
[HOE]. The models under consideration in this work are studied in the multiphoton 
regime. 

If the oscillator is highly excited, then the problem simplifies. The two-level system 
interacts with the oscillator to produce a shift as before; however, the oscillator is only 
weakly impacted by the two-level system since the two-level energy in this case amounts 
to a small fraction of the total oscillator energy. In this limit the Bloch-Siegert shift 
is accurately approximated using an adiabatic model [6], and we have obtained new 
estimates for the level splitting at the anticrossings [7| using degenerate perturbation 
theory on a rotated version of the model. In the rotated problem, multiphoton 
transitions are mediated through a complicated "perturbation" operator, but in the end 
the level-splitting is obtained approximately from a simple two-state model with only 
first-order coupling. The complicated interactions of the original spin-boson problem, 
which are responsible for energy exchange of a large number of quanta, are reasonably 
well accounted for through the lowest-order coupling in the rotated version of the 
problem. 

This approach, and the resulting conceptual simplification of the problem, is 
reasonably general and very powerful. We extended the analysis to the case of a spin- 
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one system coupled to an oscillator [8], and obtained results for the Bloch-Siegert shift 
in agreement with previous work, as well as new results for the level splitting at the 
Bloch-Siegert resonances. The accuracy was comparable to that obtained for the spin- 
boson model, and it seems clear that the approach can be applied systematically to 
higher-spin generalizations of the spin-boson model as well. We also studied a different 
generalization of the spin-boson model in which a three-level system is coupled to 
an oscillator [9J. Technical issues associated with the three-level system made the 
implementation of the rotation much more challenging; however, in the end we obtained 
good results for the level shifts, and for the level splitting at the anticrossings, as long 
as the anticrossing occured away from other strong resonances which interfere. 

In this work we turn our attention to a different generalization of the spin-boson 
model in which two different two-level systems are coupled to an oscillator. This and 
related models have been studied in connection with studies of the interaction between 
two atoms and a cavity [TUl [HI H2], and quantum entanglement [131 EE]- Our 
approach is most useful in the multiphoton regime with a highly excited oscillator. In 
this case, both two-level systems experience Bloch-Siegert shifts due to their interaction 
with the oscillator, and anticrossings occur associated with energy exchange between 
each two-level system and the oscillator. The Bloch-Siegert shift for each two-level 
system in this case is very close to what would be expected if the other two-level system 
were absent, and the level splittings at the energy exchange anticrossings are not very 
different from the single two-level version of the problem. What is new in this problem 
are anticrossings associated with excitation transfer, in which the excitation from one 
two-level system is transferred to the other. In the rotated version of this model, we find 
a spin-spin interaction term which mediates these new excitation transfer transitions. 
Our attention in this work is then focused on excitation transfer, and we find that a 
simple two-state model and degenerate perturbation theory leads to reasonably accurate 
estimates of the level splittings away from other resonances. 

The excitation transfer effect in this model is very weak. To study it in a regime in 
which energy exchange resonances do not interfere, we need to work with low oscillator 
energy (which maximizes the number of oscillator quanta needed for resonance) and 
modest oscillator excitation (since the level splitting is inversely proportional to the 
number of oscillator quanta). The reason why excitation transfer is so weak in this 
model is studied using a finite-basis approximation; in which we see that destructive 
interference occurs between the contributions from pathways involving intermediate 
states. 
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We briefly examine the possibility of reducing or eliminating the destructive 
interference, in order to make the excitation transfer effect stronger. If the model 
is modified so that the coupling between the different two-level systems is through 
conjugate oscillator operators, then some of the destructive interference is removed. If 
the model is augmented with loss terms that remove the contribution of the lower- 
energy intermediate states, then the destructive interference is completely eliminated; 
the excitation transfer effect then becomes much stronger. 
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2. Model 
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We consider the model described by the Hamiltonian 

s(i) s( 2 ) 
H = AE!-^ + AE 2 ^ + huj a' l a + U 1 (a' 1 + a) 




+ U 2 {a^ +o) 




n 



X 



(1) 



There are a pair of two-level systems, with unperturbed transition energies AE ± and 
AE 2 ; and an oscillator, with characteristic energy hu . The first two- level system 
interacts linearly with the oscillator, with a coupling strength of Ui, similarly, the second 
two-level system also interacts linearly with the oscillator, with a coupling strength of 
U 2 . The spin operators are defined in terms of the Pauli matrices according to 



where the superscript denotes which two-level system is referenced. In the multiphoton 
regime, the characteristic energy of the oscillator is much less than the unperturbed 
transition energies 




(?) 



(2) 



hu < AE 



i 



huo < AE 2 



(3) 
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As in our previous work on the spin-boson problem [16] , we find it useful to consider 
a unitary equivalent Hamiltonian 



H' = UHtf (4) 



where 



U = UiU 2 = e- lXl& y e~ iX2 *v (5) 

Since cr^ commutes with cr^ 2 \ the computation is very similar to the single two- 
level problem. This unitary transform is a straightforward generalization of one used 
previously for the spin-boson problem [T71 [18] 



3.1. Dressed Hamiltonian and "unperturbed" part 
We can write the rotated Hamiltonian as 

H' = ^0 + ^1 + ^2 + ^1 + ^2 + ^12 (6) 

where the "unperturbed" part of the rotated Hamiltonian H' is given by 

/ sW / s( 2 ) 

H = ^/\El + mly^ + ^El + m 2 y^ + huj a)a (7) 



where 



a + a' 



This is similar to the unperturbed part that we obtained previously in the case of the 
spin-boson model, where now two dressed two-level terms appear instead of one. 

In the spin-boson model, and also in other problems that we have studied, the 
"unperturbed" part ifo of the rotated Hamiltonian H' results in a good approximation 
for the oscillator and dressed transition energy of the two-level systems. It can be used 
to develop approximations for resonance conditions as we discuss in the next section. 
Because of this, we have come to think of H as describing an unperturbed version of 
the dressed problem in which no interactions occur at level crossings. Viewed in this 
way, the other terms in the rotated Hamiltonian can be thought of as perturbations. 
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3.2. Perturbations involved in energy exchange resonances 

There are now two primary perturbations V\ and V% (only a single one appears in the 
spin-boson problem since there is only one two- level system in that model), each of 
which can be described by the general formula 



V; 



hu>n 



V2UA 
, A ^ ) 



2„,2 K 



1 + 



8Ufy 
AE] 



d d 
+ 



V2Uj\ 
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dy ' dy 8Uh 2 " 



1 + 



AE] 



(9) 



where 



d 
dy 



a — a) 



(10) 



In the multiphoton regime, these terms give rise to level splittings at resonances in which 
one unit of excitation of a two-level system is exchanged for an odd number of oscillator 
quanta, as discussed in [7j. 



3.3. Potential operators 

There are two small terms W\ and W2, both described through the general formula 

2 
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AEj 
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(11) 



Since the square of the spin operator s y is proportional to the identity matrix 

2 
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u 



a 




(12) 



these terms become simple potentials 

2 
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SUfy 
AE] 



(13) 



In the large n limit which is of interest in this paper, these potentials are very small, 
and do not contribute in a significant way to either the occurrence of resonances or the 
level splitting at resonance. We neglect them in what follows. 
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3.4- Spin-spin interaction 



Finally, we find a spin-spin operator V12 given by 
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This operator has no analog in the spin-boson model. In the multiphoton regime this 
term contributes to the level splitting associated with resonances involving excitation 
transfer (where one unit of excitation in one two-level system is exchanged for one unit 
of excitation in the other two-level system, accompanied by the exchange of an even 
number of oscillator quanta). As this effect is new in this model, it will be the focus of 
our attention in this work. 
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4. Energy levels and resonance conditions 

The time-independent Schrodinger equation for H is 
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E 



flUJr 



hue 



d 2 



-— + y 2 U> + ^AE 2 + 8U 2 y 2 ^ 



(2) 



+ ^AEl + 8U 2 y 2 j-y 



(15) 



4-1. Product solutions 

We can reduce this Schrodinger equation to a purely spatial problem in y by assuming 
a product wavefunction for \1/ of the form 



* = u mum2 (y)\s,m 1 )\s,m 2 ) 



(16) 



where the \s,m,j) are spin 1/2 eigenkets (with s = 1/2 and = ±1/2). This leads to 
a one-dimensional Schrodinger equation 

d 2 



eu(y) = 



dy' 



+ v(y) 



u(y) 



where the normalized energy eigenvalue e is 
2E 

+ 1 



hu 

and where the normalized nonlinear potential v(y) is 

v(y) = ^^^ + 8^ + ^ 

We have suppressed the subscripts mi on u{y). 



AE 2 + 8U 2 y 2 



(17) 



(18) 



(19) 



4-2. Approximate solutions 

In the large n limit, the oscillator is highly excited and we would not expect the two- level 
systems to have much impact on the oscillator. In this case, a reasonable approximation 
for the oscillator is to use eigenfunctions of the simple harmonic oscillator 

u{y) = (y\n) (20) 

where \n) are the harmonic oscillator eigenfunctions. Within this approximation, we 
may write 

En, mi ,m 2 = AE 1 (g 1 )m 1 + AE 2 (g 2 )mi + nhu (21) 
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where the dressed transition energies AEi(gi) and AJ5 2 (<72) are given by 



AEjigj) = AEj ( n 



1 + 



8Ufy 2 



n) (22) 



The dimensionless coupling constants #i and #2 are defined according to 

9i = ^# (23) 

We interpret AEj(gj) as the "dressed" transition energy of the jth two-level system. 
Numerical solutions for the original Hamiltonian H are well described by this expression 
away from resonances, similar to what we found in the spin-boson problem [7J. This 
also allows us to meaningfully label our eigenfunctions u(y) with n (again suppressing 
the weak mi and m 2 dependence) 

u(y) = u n (y) (24) 
4-3. Energy exchange resonances 

Within this approximation scheme, we can develop resonance conditions for energy 
exchange resonances. For energy exchange between the first two-level system and the 
oscillator, the resonance condition is 

AEx(gi) = AnfiujQ (25) 

with An odd. These resonances correspond to the ones we analyzed using a similar 
approach in the spin-boson model [7] . The level splittings at the associated anticrossings 
are illustrated in Figure Q] under conditions where the energy exchange resonances 
dominate (open circles). Also shown are approximate results using degenerate 
perturbation theory based on the eigenfunctions of the rotated Ho problem. In this 
calculation, one state is chosen in which the first two-level system is excited, the second 
is in the ground state, and n quanta are present in the oscillator; the other state is 
chosen so that both two-level systems are in the ground state, and n + An quanta are 
in the oscillator. The level splitting at the anticrossings using degenerate perturbation 
theory is [7j 

SE min = 2|(\l/ nil /2 -l/2|Vi|\l/ n+ An,-l/2-l/2)| (26) 

One sees that the level splittings are in reasonable agreement with this approximation. 
This is similar to what we found previously in the spin-boson model. Moreover, 
the level splitting found in this example matches that for the equivalent spin-boson 
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Figure 1. Level splittings associated with energy exchange anticrossings, for a model 
with AEi = 11 huj and AE 2 = 13huj 0} n = 10 4 . Solid circles - results of direct 
numerical calculation of the original H problem; open circles - results from degenerate 
perturbation theory using eigcnfunctions of Ho- The dimcnsionless coupling coefficient 
<?2 is fixed at 0.5 for these calculations. The number of oscillator quanta exchanged 
An is indicated for each anticrossing. 



model (obtained by removing the second two-level system) such that they could not be 
distinguished if plotted together in this figure. 

Energy exchange between the second two-level system and oscillator occurs similarly 
when the resonance condition 



is satisfied. 

4-4- Excitation transfer resonances 

This system supports another kind of resonance that is not present in the spin-boson 
model. Excitation can be transferred from one two-level system to the other (along with 
energy exchange with the oscillator), which has motivated us to refer to the associated 
resonances as "excitation transfer" resonances. Consider the resonance between one 
state, with energy 




(27) 



E 



'no -1/2,1/2 





and another, with energy 



E no+ An, 1/2,-1/2 



-AE( 9l ) - -AE(g 2 ) + hu (n + An) 



(29) 
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Figure 2. AE2 = l5Tiu>o, tiq = 9600, and An ranging between -10 and 10. The 
lines associated with An = — 10 and An = 10 are denoted, with the lines in between 
corresponding to even An. 



We obtain the resonance condition by requiring the two basis states to have the same 
energy 

Eno, -1/2,1/2 = -Eno+An.,1/2,-1/2 

(30) 

This is consistent with the constraint 

AE(g 2 ) - AE(g 1 ) = An?ko (31) 

with An even. 

Results from a computation based on the WKB approximation of resonance 
conditions are illustrated in Figure [2j The WKB approximation for this case is discussed 
in Appendix A. One observes that as long as AE 2 is greater than AEi + AnhuiQ, for 
every g 2 a value can be found for gi at which a resonance occurs. Similarly, as long as 
AE\ is greater than AE 2 — Anhu , for every g 1 there occurs a value of g 2 at which a 
resonance occurs. 
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As we observed in the case of a single two-level system, energy levels split at anticrossing 
when resonances occur. We can compute this splitting in at least two ways: either 
by direct numerical diagonalization of the full unrotated Hamiltonian H, or by using 
degenerate perturbation theory on the rotated Hamiltonian. 

5.1. Degenerate perturbation theory 

In the vicinity of a level anticrossing that is free of energy exchange resonance disruption, 
we can investigate the level splittings using a two-state approximation of the form 

V> = c o0o + ci0i (32) 

0o = Mn(y)|s,-l/2)|s,l/2) 0i = Un+ A„(y)|M/2)|s,-l/2) (33) 

The two-state problem leads to the following characteristic equation for the energy levels 

E ( C » ) = ( £ » <*>'^>> ) (34) 

Since we are at resonance 
E = Ei 

and the level-splitting is given by given by 

SE mm = 2|(0 o |Vi 2 |0i)| (35) 

5.2. Level splitting estimates 

The matrix element that appears here can be written as 

(0o I V u 1 0i ) = 2 ^o^^^^A2 = In (36) 



where the integral I12 is given by 



'12 



u n (y) 



1 



1 + 



1 + 2y 



AEI 



u n+ An(y)dy (37) 



AEf 

We can compute this integral using either numerical wavefunctions obtained by 
diagonalizing H , or by utilizing the WKB approximation [T5] . 

The level splitting on resonance can be expressed in terms of the integral J12 as 

5E min = °9i92 (38) 

n 
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9i 

Figure 3. Level splittings for AEi = 53 Uujq and AE2 = 51 huio, no = 1000, and 
An = —2. Solid circles - results of direct numerical calculation of the original H 
problem; solid line - results from degenerate perturbation theory using eigenfunctions 
of Hq. The results are shown as a function of gi, where 52 has been optimized to 
minimize the level splitting. 

The results of the level splittings from a direct numerical computation using the 
unrotated H and degenerate perturbation theory (using the WKB approximation) are 
illustrated in Figure [3J We can see that the perturbation theory results are in excellent 
agreement with the exact result. In this example, we have selected large two- level system 
transition energies relative to the oscillator energy in order to minimize the impact of 
energy exchange resonances. In addition, we have chosen a moderate value for n (instead 
of a larger value) to increase the level splitting. 
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The level splitting associated with an excitation transfer resonance is a weak effect in 
this model. We can better understand the slow dynamics of the excitation transfer by 
making use of a finite basis expansion. One can see from this kind of calculation that 
destructive interference between different pathways produces a very small second-order 
coupling between initial and final states. If this destructive interference can be broken, 
then the second-order coupling is greatly increased. 

6.1. Finite-basis approximation 

Indirect coupling between initial and final states associated with an excitation transfer 
process can be analyzed in the weak coupling limit through the use of a finite-basis 
approximation. Consider a finite-basis approximation with six basis states 

* = E Cj % (39) 

where the basis states are 

$i = |n)|a,l/2)|a,-l/2) 
*2 = |n-l)|a,-l/2)|s,-l/2) 
*3 = |n+l)|a,-l/2)|s,-l/2) 
*4 = |n-l>|M/2>hl/2> 
$ 5 = |n + l}|s,l/2)|s,l/2) 

* 6 = |n)| a) -l/2)|a,l/2) (40) 



The energy levels and coupling are indicated schematically in Figure HI The excitation 
transfer process in this case would take the system from an initial state $i (with an 
excited first two-level system, and a ground state second two-level system) to a final 
state $6 (with a ground state first two-level system, and an excited second two-level 
system). Since the Hamiltonian H does not couple the basis state $i to the basis state 
$6 directly, the coupling between the two states is indirect. The intermediate states 
$2, • • • $5 will provide the dominant pathways between the initial and final states in the 
case of weak coupling. 

6.2. Indirect coupling 

It is possible to obtain an indirect interaction between basis state $i and basis state 
$6 by eliminating the intermediate states algebraically, as we illustrate in what follows. 



Excitation transfer in two two-level systems coupled to an oscillator 



16 





Figure 4. Schematic of levels for a pair of two-level systems with indirect coupling 
between two degenerate states. The basis state index j is indicated near the energy 
level. Boxes are included that give the number of oscillator quanta; arrows indicate 
whether the first two- level system is excited or not (first arrow), and whether the 
second two- level system is excited or not (second arrow). 



The finite-basis equations for the expansion coefficients c\ ■ ■ ■ c 6 can be written as 



Ec\ = H\C\ + U\\fnc2 + U\Vn + IC3 + L^V^Q. + t^Vn + IC5 

Ec 2 = H 2 c 2 + U\\fnc\ + U2^/ncQ 

Ec 3 = H 3 c 3 + UiVn + lei + U 2 Vn + lc 6 

= H4C4 + U 2 Vnci + U\ \pnc<o 
Ec 5 = F5C5 + U 2 Vn + lei + UWn + lc 6 

Ec% = Hqcq + \J2\fnc2 + U 2 Vn + IC3 + U\y/ncA + f/i\/n + IC5 (41) 



The algebraic elimination of the coefficients C2 through C5 produces two coupled 
equations of the form 

Ec 1 = [H 1 + Y.i(E)]c 1 + V 1 s< < E)c % 

Ec 6 = [H 6 + S 6 ( J B)]c 6 + V 61 (E)a (42) 



where 



V 16 (E) = V 61 (E) = U1U2 



n 



n+1 
+ ^ + 



n 



+ 



n+1 



(43) 



— H 2 E — H 3 E — H4 E — 
The self-energy terms ^i(E) and T*g(E) are small in the case of weak coupling, and 
contribute to the energy at which the resonance occurs, but not to the splitting. 
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6.3. Level splitting 

Within this model, the resonance condition is 

Hr + H^E) = H 6 + Z 6 (E) (44) 

which is the finite-basis approximation equivalent to Equation ( 13TI) when no oscillator 
quanta are exchanged in association with the excitation transfer process. The level 
splitting at resonance in this finite basis model is 



5E mm = 2yjV 16 (E)V 61 (E) (45) 

If the self-energy terms can be neglected, then the energy eigenvalue E is very nearly 
equal to the unperturbed level energy of the initial and final states 

E — Hx — H 6 (46) 



V 16 (E) = ,°^ 2 (47) 

ibV ' / A N2 ff. s 2 V 1 



We can evaluate the indirect coupling term Vi 6 (E) at resonance to be 

2hu U 1 U 2 
(AE,) 2 - (huj o y 

where AEi must approximately the same as AE 2 for the resonance condition of Equation 
( I46p to be satisfied. When the characteristic oscillator energy tkvo is much smaller 
than the transition energy AE\, then the level splitting on resonance 5E m i n in this 
approximation evaluates to 

SEmin = ^Mlii m 

n 

This is equivalent to the splitting we obtained using degenerate perturbation theory in 
the rotated frame [Equation ( 1381) ] as long as the integral I\ 2 is unity. From an inspection 
of Equation (1A.9I) . we see that this is the case as long as 

9i,92 < 1 (49) 

Our result obtained in this section then is the weak coupling limit of what we obtained 
above. The generalization to a multi-mode system is considered in Appendix B. 

6.4- Destructive interference and loss 

The second-order coupling coefficient Viq(E) for indirect coupling between state <J>i and 
$6 is much smaller than the direct coupling coefficients with intermediate states. For 
example, we may write 

Vw(E) 2g 2 huj 

U x Jh ~ n AEt 1 j 
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This reduction of coupling strength is due to destructive interference between the 
different contributions in Equation ( H3l) that make up Viq(E). 

To show that this is so, we consider a modified version of the model in which the 
destructive interference is removed. Consider the two-spin plus oscillator Hamiltonian 
augmented with a loss term —ihf(E)/2 

sW , ih~ 

H l = AE 1 ^- + AE 2 ^- + hu ffa--T{E) 

+ Z7i (a f + a) —2- + U 2 (a f + a) (51) 

This loss term accounts for energetic decay processes of the low-lying intermediate states 
at the transition energy of the two-level systems. Only intermediate states with both 
two-level systems in the ground state ($2 and $3) can decay in this model, since these 
states have unperturbed energies much less than the available energy E. The system 
cannot decay in such a way as to produce $4 and $5 as final states. 
The coefficients c 2 and C3 in this kind of model now satisfy 

H 2 — \c 2 + U\\fnc\ + U 2 ^fncQ 

( ihT\ 

Ec 3 = lH 3 - — \c 3 + UWn + lei + U 2 Vn + lc 6 (52) 

The indirect coupling coefficient Vie(E) in this model is changed due to loss effects. We 
may write 



V 16 (E) = XJJJ 2 



n n + 1 n n + 1 

+ - „ , — 7T + — — + 



(53) 



E - H 2 + ihT/2 E-H 3 + ihT/2 E - H 4 E - H< 
In the limit that the loss term becomes very large then we may write 

2U 1 U 2 AE 1 n 

Vl ^ E) = ~ (A^) 2 - (^o) 2 ^ " °° ) (54) 

The ratio of this indirect coupling coefficient (which is now free of destructive interference 

effects) to the direct coupling coefficient between states $1 and $2 becomes 

m - "2 92 (55) 

The removal of destructive interference increases the indirect coupling coefficient by a 
factor of 

|Vi 6 (£)|r=oo = nAEi 
\Viq{E)\t=q huj 
In Appendix C, we consider a different modification of the model in which the 

coupling between the different two- level systems and oscillator is through conjugate 

oscillator operators. Part of the interference is removed in this model. 
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7. Summary and conclusions 

In previous publications, we studied the shifts and splitting of energy levels in the 
multiphoton regime in the spin-boson model, and in generalizations of the spin-boson 
model to the case of an oscillator coupled to a spin-one system, and to a three-level 
system. In these works, we made use of a unitary transformation which rotates the 
Hamiltonian into a form in which terms primarily responsible for energy level shifts are 
separated from terms responsible for the level splitting. Here the same general approach 
is applied to a model involving a pair of two-level systems coupled to an oscillator, where 
the transition energies and coupling strengths can be different. The energy levels can be 
approximated accurately using a WKB approximation based on the unperturbed part 
of the rotated Hamiltonian H , similar to what we found previously. The level splitting 
at the Bloch-Siegert anticrossings are also described accurately (away from anticrossing 
resonances) using degenerate perturbation theory based on the V\ and V2; the situation 
is very similar to what we found in the spin-boson problem. 

What is new in this model (with no analog in the spin-boson problem) is the 
excitation transfer effect, in which a transition in one two-level system occurs in concert 
with a transition in the other two-level system. Approximate level energies derived from 
the unperturbed Hamiltonian H can be used to locate excitation transfer resonances 
(see Figure [2]). The excitation transfer effect is mediated by the spin-spin term V12 in 
the rotated Hamiltonian; level splittings away from the Bloch-Siegert (energy exchange) 
resonances are accurately modeled using degenerate perturbation theory based on the 
spin-spin term. The effect is weak in this model due to destructive interference effects. 
Destructive interference was examined in the context of a finite-basis approximation 
appropriate to the weak coupling limit of the model, in which contributions from 
different paths can be seen to cancel. A version of the model augmented with loss 
in the intermediate states removes the destructive interference, and leads to a drastic 
increase in the indirect coupling term. 

A consequence of the destructive interference is that the excitation transfer effect 
is independent of mode excitation in the weak coupling limit. If there are many modes 
present, then it is likely that interference between the contributions of the different 
modes will further decrease the effect, especially if the atoms or molecules represented 
by the two-level systems are far apart. This effect is discussed in Appendix B. Hence, 
one would expect excitation transfer to be a weak and short range effect in a multi-mode 
system where the destructive interference involving different pathways is not removed. 
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Such is the case in phonon-mediated excitation transfer. 

One question raised from the analysis and discussion presented here concerns 
the possibility of developing an excitation transfer system in which the destructive 
interference is reduced or eliminated. In this paper we have noted a reduction of 
interference in a modified version of the model with conjugate coupling (Appendix 
C), and elimination in a model with large loss in intermediate states (Section 6). In 
either case, the excitation transfer rates and level splitting will be greatly increased. A 
conjugate coupling scheme could be developed in an electromagnetic resonator under 
conditions where one two-level system couples to the electric field, and the other two- 
level system couples to the magnetic field (since electric and magnetic fields in a 
resonator are conjugate variables). To implement a physical system in which loss impacts 
the destructive interference, one would require a strongly driven mode which sees low 
loss at the resonant frequency, but is very lossy at higher frequency corresponding to 
the transition energies. 
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We have found the WKB approximation to be very effective for the large n limit of this 
problem. Within the WKB approximation, the eigenfunctions u(y) are written in the 
form [19] 

. , ^sin0 

u(y) = C—^ A.l 



where 

d 



, 0(2/) = v(y) (A.2) 
dy 



V (y) = ^e-v(y) (A.3) 

and where C is a normalization constant. Written in terms of WKB eigenfunctions, the 
integral 1 12 becomes 

T - r C f° sin[0o(y)]sin[^i(2/)] 



, 8C/ 2 V 



«W 1 + ^ 1 + 



A£ 2 y y ae| y 

The product of sine functions in the numerator can be decomposed into terms 
which involve rapid oscillations, and slow oscillations. For this, we can make use of the 
trigonometric identity 

• (X \ ■ (A\ COS(01 - O ) - COS(01 + O ) , 

sm(0o) sm(0i) = - (A.5) 

On the RHS, the first cosine function has a slowly varying phase, and the second one 
has a rapidly varying phase (especially if n is very large). Usually the contribution from 
the cosine with the slowly varying phase dominates, in which case we may write 

T OA roc co S [A^)] , (A R . 



AE 2 J \ AE 2 J 

In writing this we assume that the two WKB momentum variables rjo and rji are not 
very different 

vo(y) ~ vi(y) -> v(y) (A-* 7 ) 

The normalization constants Co and C\ in this approximation satisfy 
C 2 r 00 1 

dy = 1 (A.8) 



2 J-00 rj(y) 
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Table Al. Results for I12 at resonance for AEi = 11 Tvjjq and AE2 = 15 Tiloq for 
no = 9600. Values of the dimensionless coupling constants g\ and (72 at resonance (first 
two columns); integral computed from eigenfunctions of Hq (third column); integral 
computed from the WKB approximation (fourth column). 



9i 


92 


\h2\ 


\I l2 [WKB]\ 


0.50 


0.1611206 


0.2224 


0.2235 


0.60 


0.2734154 


0.1946 


0.1951 


0.70 


0.3666892 


0.1733 


0.1734 


0.80 


0.4528576 


0.1562 


0.1562 


0.90 


0.5353554 


0.1421 


0.1420 


1.00 


0.6156505 


0.1303 


0.1301 



This allows us to recast the WKB approximation in the form 

cos[A0(y)] 



-dy 



AE? I \ AE\ ) 

'■vi ,00 / - 1 — (A.9) 

dy 



3 v(y) 

In Table IA1I we give the results of computations of the magnitude of the integral 
I/12I from the numerical solution of Equation (IT7|) . and from the WKB approximation 
of Equation (IA.9I) . One sees that the WKB approximation in this case is very close to 
the numerically exact result. 
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There is an additional mechanism that produces destructive interference in the multi- 
mode generalization of this problem. In particular, interference can occur between the 
contributions of the different modes. 

Appendix B.l. Multi-mode Hamiltonian 

For simplicity, suppose that we have many modes that can be described in terms of 
wavevectors k; and that the mode frequencies are a function of the wavevectors according 
to the dispersion relation 

L) = u>(k) — > u; k (B-l) 

The multi-mode generalization of the two-spin problem then might be described 
according to a Hamiltonian H mm of the form 



„(1) .(2) _^ 

AE 1 ^- + AE 2 ^- + Y, ^ k a k a k + £ 



n 



C/i(k)a k e 4k - ri + C/ 1 *(k)aj c e 



t „— ik-ri 



2# 



+ ]T C/ 2 (k)a k e ik - r2 + C/ 2 *(k)a k e 



t „-ikr 2 



2j? 

h 



(B.2) 



where ri and r 2 are the position vectors of the two atoms. 



Appendix B.2. Indirect coupling without loss 

It is possible to develop a finite-basis model similar to that discussed in Section 6, but 
now for the multi-mode case, to obtain the indirect coupling term for the weak-coupling 
limit. If we carry out such a calculation, in place of Equation (1471) . we obtain 



V ie (E) = 

(AE,) 2 - (hcoo) 



L/i(k)05(k)e* , ( ri - ra ) , 3 , f C/ 1 *(k)£/ 2 (k)e- ik < ri - r2 ) , 



d k - / 1V \ " ' r rf*k (B.3) 

AE\ - huj k J AE ± + huj k 

The contribution of a single highly off-resonant mode appears on the same footing as 
contributions from many other off-resonant modes, each weighted by a phase factor. 
In the event that the separation |ri — r 2 | is large, then severe cancelation can occur 
between the contributions of the different modes; hence we would expect interactions 



Excitation transfer in two two-level systems coupled to an oscillator 24 

to be of short range in the multi-mode case. Indirect coupling through coupling to 
common phonon modes would be described using such an approach (see for example 
[20] . If one adopts for the multi-mode oscillator longitudinal photon modes (all with 
zero energy), then indirect coupling between electric dipoles can be thought of in this 
way; with an overall |ri — r2 1 — 3 dependence of the dipole-dipole matrix element in the 
absence of retardation for single longitudinal photon exchange. Excitation transfer in 
this case is known as resonance energy transfer in biophysics |21j . 

In general, any long range interactions in this kind of model will be dominated by 
resonant interactions if resonant modes exist. Otherwise, the indirect interaction will 
most likely occur only at short range. 

Appendix B.3. Indirect coupling with loss 

If some mechanism is present that can remove the destructive interference effect in 
the different pathways, then it may be possible to reduce or eliminate the destructive 
interference effect in the multi-mode version of the problem. The Hamiltonian in this 
case can be written as 



-(2) hf(E) 

H mml = AE l ^+^E 2 ^-+Y,^\A.^- i —ir-+Y J 
ft ft k i k 



!7 1 (k)a k e lk - ri +C/ 1 *(k)a k e 



t „— ik-ri 



2§x 



^ 2 (k)a k e 4k - r2 + C/ 2 *(k)a k e 



ftt „— ikr 2 



2^ 
ft 



(B.4) 



An analogous finite-basis type model can be developed for this kind of model, leading 
to an indirect coupling in the r — > oo limit of 



V 16 (E) = ™ lU2 AE in 



(AE x f - (ft^o) 2 

t/ 1 (k)C/ 2 *(k)n(k)^ k ^-^) d 3 k r UZ(k)U 2 (k)[n(k) + l]e- ik '' r '- r ^ 3 

AEi - ftu; k J AEi + ftw k { ' ' 

Things are qualitatively different in this case, since the indirect coupling now depends 

on the number of quanta n(k) in the different modes. If a single off-resonant mode is 

very highly excited, then the interaction can take place over a much longer range since 

the unexcited mode can no longer destructively interfere. 
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It is possible to eliminate some of the destructive interference by allowing the two 

different two-level systems to couple with the oscillator through interactions that do not 

commute. Consider a modified version of the model described by the Hamiltonian 

s(2) 2 ~(i) , s (2) 

H e = AE 1 ^ + AE 2 ^ + hu tfa + U 1 (tf + a)^ + iU 2 {tf-a)—?-{C.l) 

If adopt the same finite-basis approximation for the wavefunction \1/ as discussed in 
Section 6 

* = £ c i *i (C2) 
i=i 

then the expansion coefficients satisfy 



Ec\ = H\c\ + U\^/nc 2 + Uxvn + IC3 + iU 2 y/nc4 - iU 2 \fn + IC5 

-Ec 2 = + Uiy/nci + iU 2 Vnc e 

Ec 3 = H 3 c 3 + £/"iVn + lei - iU 2 y/n + lc 6 

SC4 = #404 — iU 2 ^fnc\ + Uiy/ncQ 

Ec 5 = H 5 c 5 + iU 2 y/n + lc\ + Ux\/n + Icq 

Ecq = HqCq + iU 2 ^fnc 2 - iU 2 y/n + lc 3 + L r i V / ^C4 + f/iVn + lc 5 (C.3) 



i?9 



+ 



The indirect coupling coefficient Vie(-E) in this case is 

V 16 (E) = zU.U, 1 U H + 1 

At resonance, we assume that 
E = H\ = H§ 

to obtain 



■n 



n + 1 



2iU 1 U 2 hujQ{2n + : 



(C.4) 



(C.5) 



(C.6) 



This indirect coupling is much larger than we obtained with the original model, and 
also depends on the mode excitation. However, we see that it is smaller in magnitude 
by Hujq/AEi than the result [Equation (1541) ] we obtained in which the destructive 
interference was completely eliminated by loss. Hence, we have removed some of the 
destructive interference with conjugate coupling. 
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